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' We study the fluctuations that are predicted in the autocorrelation 

function of an energy eigenstate of a chaotic, two-dimensional billiard 
by the conjecture (due to Berry) that the eigenfunction is a gaussian 
random variable. We find an explicit formula for the root-mean-square 
O \ amplitude of the expected fluctuations in the autocorrelation function. 

These fluctuations turn out to be 0(h 1 ^ 2 ) in the small h (high energy) 
limit. For comparison, any corrections due to scars from isolated periodic 
orbits would also be 0{% 1 ^ 2 ). The fluctuations take on a particularly 
^ \ simple form if the autocorrelation function is averaged over the direction 

of the separation vector. We compare our various predictions with recent 
numerical computations of Li and Robnik for the Robnik billiard, and 
find good agreement. We indicate how our results generalize to higher 
dimensions. 
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Two-dimensional billiards which are classically chaotic have proven to be an efficient 
laboratory for the study of quantum chaos. The energy eigenvalues and eigenfunctions 
can be computed with good accuracy, and compared with theoretical predictions of their 
properties. These predictions are typically semiclassical in nature, involving properties that 
are expected to be emergent in the formal limit of h — *■ 0. 

In practice, numerical methods are used to find the eigenvalues k 2 and eigenfunctions 
"0fc(x) of the time-independent Schrodinger equation, 

(V 2 + fc 2 )^(x) = 0, (1) 

where x is in the domain B of the billiard, and the Dirichlet boundary condition 

^ fc (x) =0 on dB (2) 

is imposed. Then large k corresponds to small h; an expansion of some quantity in powers 
of the wavelength A = 2n/k corresponds to an expansion in powers of h. 

Our focus here will be on the autocorrelation function C^ ^s), introduced by Berry 
0]. Given an eigenfunction ^(x), a separation vector s, and an averaging region R, the 
autocorrelation function is defined to be 

Ab f 

C k ,n{s) = — / d 2 x ^fc(x + |s)V>fe(x - |s) , (3) 

A R JR 

where Ab is the area of the billiard and Ar is the area of the averaging region R. The 
eigenfunction ^fc(x) is assumed to be real, and normalized in the usual way, 

/ d 2 x ^ fc 2 (x) = 1 . (4) 

However, our normalization of Cfc^s) differs slightly from Berry's; we will discuss the reason 
for this later. 

Berry conjectured that an energy eigenfunction in a chaotic billiard would appear locally 
to be a superposition of plane waves with random directions of the momenta and random 
phases. This implies that the expected value of C^^s) is 

Jd 2 p6(p 2 -k 2 )e^ s 

{CkAs)) ~ TW^T 

= Mks) , (5) 

where s = |s|, and Jo(x) is a Bessel function. We now know that there are corrections to this 
result ("scars" @) associated with isolated periodic orbits in the classical billiard. Assuming 
that the averaging region R encompasses many wavelengths in the direction perpendicular 
to each orbit giving a scar, then the scar corrections to @ are suppressed by Oifi 1 ^ 2 ). 

In the limit that the number of superposed plane waves becomes infinite, the central 
limit theorem tells us that the function ^(x) can be treated as a gaussian random variable 
|T|J^Q]. This means that we have prior information (in the sense used in Bayesian statistical 
analysis) about ^(x) which can be represented by a functional probability distribution of 
the form 
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P(V>*) = ATexp 



-- / d 2 x 1 f d 2 x 2 V , fc(xi)fr(x 1 ,x 2 )^fe(x 2 ) 
2 Jb Jb 



(6) 



Here A/" is a normalization constant, and X(xi,x 2 ) is the inverse of Ag 1 Jo(/c|xi — x 2 |) in the 
sense that 



/ gPx 2 ^(xi, x 2 ) J (fc|x 2 - x 3 |) = A B 5 2 {^i - x 3 ) 
Jb 



(7) 



with X! and x 3 restricted to B. The angle brackets in @ are then defined to represent an 
average over the probability distribution (|6|). Thus we have 



(^ fc (x 1 )'0 fc (x 2 )) = J (fc|xi - x 2 |) 



(8) 



Combining this with the definition (|3|) of C kjR (s) gives (C^^s)) = J (ks), in agreement 
with (0). However, the probability distribution (Q) also contains information about the 
fluctuations of Cfc^s) about (Ck }R (s)). Our goal is to study the properties of these 
fluctuations. 

Before proceeding, let us recall that there is striking numerical evidence in favor of 
another consequence of (^]): specifically, the probability P(x)dx that Vfc( x ) has a value 
between \ an d X + a t a particular point x is given by 



P( X ) = (A B /2vr) 1 / 2 exp(-|W 



(9) 



This prediction can be tested by dividing the billiard into small pixels, and making 
a histogram of the value of the eigenfunction at each pixel. This was first done by 
MacDonald and Kaufman [|| in their study of eigenfunctions of the stadium billiard with 
k 2 As — 1.3 x 10 4 . More recently, Li and Robnik || studied eigenfunctions of the Robnik 
billiard || with )z 2 Ab — 2.5 x 10 6 , and found excellent agreement with (^). Generally, the 
prediction (|Bp is expected to be valid provided that distortions of the billiard boundary on 
the scale of the wavelength A = 2njk do not permit the formation of an integrable billiard 

Our main tool in studying the fluctuations of C k r(s) about (C k _r(s)) will be the relation 

(^fe(xi)^fc(x 2 )^&(x3)V'fc(x4)) = (^ fc (xi)^fc(x 2 ))(^ fc (x 3 )?/; A .(x4)) 

+ (^fc(xi)^ fc (x 3 ))(V' fc (x 2 )V' fc (x4)) 

+ (^ fc (xi)^(x 4 ))(^ fc (x 2 )^ fc (x 3 )) . 



(10) 



In particular, we consider the quantity 

Afc,R( s l5 S 2) = ( Cfc p(Si) 



(C k , R (si)C k , R (s 2 )) - (C k , R { Sl ))(C h , R {s 2 )) 



(11) 



A kjR (s, s) ' represents the root-mean-square discrepancy to be expected between Ck >R (s) 
and (Ck >R (s)) J7|, while A fc J? (s 1 ,s 2 ) tells us whether the discrepancies for s = Si are 
correlated with those for s = s 2 , and whether this correlation is positive or negative. 
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Let us note that quantities such as ( | V'fc ( x i ) 1 2ri I V'fc ( x 2 ) | 2t " ) have been computed previously, 
but with the angle brackets representing an average over a random potential [jj. This 
random-potential average was subsequently shown to be equivalent to the average over the 
eigenfunction probability distribution P(ipk) UTOfl - 

Returning to ([□]), we use the definition @ of Ck,ii{s) and the combinatoric property 
(0) to get 

A 2 ( ( 

A fc ,fl(si, s 2 ) = -t# / d 2 xi / d 2 x 2 C0fc(xi + |si)V>fc(x 2 + §s 2 ))(Vfc(--i - isi)^ fc (x 2 - |s 2 )) 



A R Jr 



R 



+ + |si)^ fc (x 2 - is 2 ))('0fc(x 1 - isi)^ fc (x 2 + |s 2 )) 



;i2) 



Now using (H), we find 



A fc ,fl(si,s 2 ) = -to- / d 2 £i / d 2 x 2 J (&|u + s_|) J (&|u - s_|) 



(13) 



(14) 



+ J (k\u + s+|) «7 (fc|u - s+|) 

where we have defined 

u = xi — x 2 and s± = |(si ± s 2 ) . 

To proceed further, we assume that the area An is large, in the sense that both 

Ar > A 2 and A R > s 2 2 . (15) 

In this case, the argument of each Bessel function is large over most of the range of the 
integrand, and we can use the asymptotic formula 



J (x) 



2 \ V2 / n \ 

COS X — 7 

71X I v 



(16) 



which in fact is an excellent approximation for all x > 1. Making the replacement (|i~6"D, 
expanding in s/u, and keeping only those terms which are not suppressed by extra powers 
of either ku or s 2 /u 2 , we have 



Jo(k\u + s|) Jq(A;|u — s|) ~ — — cos(ku + ks cos 6 — |J cos(ku — ks cos6 — | 

sin(2A;M) + cos(2/cscos6 



irku 
1 



(17) 



where 6 is the angle between u and s. We now use ( |TTD in (|13"D, and notice that the rapid 
oscillations of sin(2/ru) will cause this term to integrate to zero (to a good approximation). 
Thus we find 



A 



fc,H(Sl,S 2 ) 



7rkA 2 R Jr 



d 2 x± / d 2 X2U 1 cos(2/cs_ cos^_) + cos(2/cs + cos^ 



R 



(18) 
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where 9 + is the angle between u and s + (s_). 

To get a more explicit formula, we need to choose the shape of the averaging region R. 
For a disk of diameter d and area A R = ^nd 2 , the integrals in ( |T8| ) can be done in closed 
form by changing the integration variables to u = Xi — x 2 and v = Xi + x 2 , integrating over 
v subject to the constraints |v ± u| < d, and then integrating over the magnitude of u to 
get 



A fc , fi (si,s 2 ) 
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^ d9 



^/ 2 kA% 2 Jo 2vr 



cos(2A;s_ cos 9) + cos(2/cs + cos i 



(19) 



where 9 + and 9_ have each been shifted and renamed 9. Performing the integral over 9 gives 
us our central result, 



A fcijR (si, s 2 ) 
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3n 3 / 2 kA]( 2 



J (fc|si - s 2 |) + J (A;|si + s 2 |) 



(20) 



We now turn to a study of the implications of fl2"0"|). The expected discrepancy between 
C fcijR (s) and {C k)R (s)) is given by 



A fcifl (s,s) 1 / 2 = 1.38 {k 2 A R )-^ 



\ + \U2ks) 



1/2 



(21) 



where the function in square brackets attains its maximum value of one when ks = 0. 
Since fl2~I ) is proportional to A; -1 / 2 , it is 0{h 1 ^ 2 ); thus, the RMS amplitude of the expected 
fluctuations in Ck sR (s) vanishes in the classical limit. However, this amplitude is not 
numerically small unless A R ^> A 2 . Both of these points are in accord with Berry's original 
(qualitative) discussion of the approach of Ck r(s) to (Ck R (s)) as h — > Furthermore, 
A fcii? (s,s) 1 / 2 is the same order in h as any corrections due to scars. This is consistent with 
the idea H that scars represent a particular organization of the gaussian fluctuations in the 
eigenfunction, rather than constituting an additional phenomenon. 

For comparison, we turn to the numerical results of Li and Robnik || for the Robnik 
billiard 0. We computed Ck, R (s) using the eigenfunction with k = 790.644, shown in fig. (3) 
of 0, which was kindly supplied to us by Li and Robnik. The averaging region R was taken 
to be a disk with diameter d = 0.273 = 34.4 A. For this value of d, the coefficient of the 
bracketed function in ([H]) is 0.100. The leading corrections to ( p0|) from terms we have 
neglected (due to our various approximations) are suppressed by an extra factor of either 
1/kA 1 ^ = 0.005 or s 2 /A R ; for ks = 30, s 2 /A R = 0.025. In figs. (1-12), we plot the actual 
correlation function Ck ;R (s) as a solid line, along with a shaded band encompassing the range 
(Ck, R {s)} ± Ak, R (s, s) 1 / 2 . The inset in each figure depicts the Robnik billiard, and the filled 
circle shows the averaging region which was used. The double-headed arrow in each inset 
has unit length, and its direction shows the direction of s. We see that the actual Ck >R (s) 
always lies within the shaded band for a majority of the time, but does have (sometimes 
large) excursions outside of it. Without attempting a detailed quantitative analysis, we can 
say that these graphs are qualitatively consistent with what we expect. 

Li and Robnik |J suggested that the discrepancy between Ck, R (s) and (Ck jR (s)) could 
be reduced by averaging C^^s) over the direction of s. Let us define 



C, 



k,R{S 



2tt 



g C*,*W)) 



(22) 
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where s(0) = (scos0, ssin0). Obviously, we have 

(C k , R (s)) = (C fc> ji(s)> = Jo(ks) . (23) 

We also define 

A k}R (si,s 2 ) = (C k}R (si)C ktR (s 2 )) - (C kjR (si))(C ktR (s 2 )) . 
Jo 2tt Jo 2tt 



Then, using (20), we find 

32 r 27T dS 

A k , R ( Sl ,s 2 )= / f J (k[sl + si ±2s lS2 cos ^ , (25) 

where <fi = <\>\ — <p 2 ; this integral can be performed to yield 

32 

A fc ,fl(si, s 2 ) = — — ^ J (/csi) J (/cs 2 ) . (26) 



3tt 3 / 2 A;A 



R 



The fact that A fc j ? (si,s 2 ) is proportional to (C k>R (si)) (C k>R (s 2 )) has dramatic conse- 
quences; it implies that C k)R (s)/ {C k>R (s)) must be independent of s. To demonstrate this, 
we choose a set of orthonormal basis functions f n (s), n = 0, 1, 2, ... , with fo(s) chosen to 
be equal to (C k , R (s)). We require orthonormality in the sense that 

rod 

/ dsw(s)f n (s)f m (s) = 5 nm , (27) 
Jo 

where w(s) is any weight function which ensures the convergence and correct normalization 
of the integral when n = m = 0. (Since (C kjR (s)) = Jo(ks), we could construct such a 
set of basis functions by starting with the Bessel functions J n {ks) and then performing 
Gram-Schmidt orthogonalization.) Once we have the basis functions, we can write 

oo 

C k , R ( s ) = £ Cnfn(s) , 
n=0 

roo _ 

c n = dsw(s)f n (s)C kR (s) , (28) 
Jo 

where the c n 's should be regarded as random variables. By construction, we have 

(c n ) = 5 n0 . (29) 
Using (|26|), we can also compute the expected value of (? n . We find 

(c 2 n ) - (c n ) 2 = ds 1 w(s 1 )f n (si) ds 2 w(s 2 )f n (s 2 )A kiR (s u s 2 ) 
Jo Jo 

32 



3n 3 / 2 kA R /2 Jo 



ds 1 w(s 1 )f n (s 1 )fo(s 1 ) / ds 2 w(s 2 )f n (s 2 )fo(s 2 ) 
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(30) 



Thus (c 2 ) = if n 7^ 0, indicating that the probability distribution has no support for any 
nonzero c n other than c . Therefore C^r(s) oc fo(s) = (Ci. } r(s)). 

However, we must remember that there are additional contributions to A^r^si, s 2 ) which 
are suppressed by an extra factor of either l/kA]( 2 or s 2 /Ar, and that these will make small 
corrections to the functional form of Afc^si, s 2 ). This means that Cj-,r(s) / (Cf-,R(s)) should 
be independent of s up to corrections of order l/kA 1 ^ 2 and s 2 /Ar. 

The discrepancy between Ck : R(s) and (Ck : R(s)) is governed by A kjR (s, s) 1//2 . In figs. (13- 
16), we plot the actual direction-averaged correlation function C*a;,r(s) as a solid line, using 
the same four averaging regions as before. We also plot a shaded band encompassing 
the range (Cfc^s)) ± Afc^s, s) 1//2 . We see that the actual Ck,R(s) is consistent with our 
expectations. In fig. (17), we plot [Ck,R(s) — Ck, R (0)Jo(ks)} + Ck, R (0) for the four averaging 
regions; this quantity should be independent of s and equal to Cfc^(0). (We plot it instead of 
the ratio C^r(s) / Jo(ks) because the latter is dominated by numerical errors near the zeros 
of its denominator.) The plots are remarkably flat; the small glitches which are present are 
most likely due to the build-up of round-off errors in the numerical computation. These 
plots confirm our prediction that Cfc,i?( s ) / (C^^s)) should be independent of s. This result 
is an incisive test of the validity of (|26|), and by implication, (|10|). 

This concludes our analysis of the fluctuations in the autocorrelation function for the 
case of a circular averaging region. We now consider the dependence of Afc^Si^) on the 
shape of the averaging region R. 

For noncircular R, the integrals in (|18|) cannot be done in closed form for nonzero s 12 . 
For the special case Si = s 2 = 0, however, we can evaluate (|18|) for a rectangular averaging 
region. This will give us the shape dependence of A k R (0,0), and therefore (we hope) some 
idea of the shape dependence of Ak, R (si, s 2 ) for general Si and s 2 . For a rectangle with edge 
lengths a and b, each of which is much greater than the wavelength A = 2ir/k, we find 

AjUO, 0) = — i_ [g(0 + giC 1 ) - (£ + T 1 ) 372 ] , (31) 
3nkA R 

where A R = ab is the rectangle's area, £ = b/a is the ratio of edge lengths, and 

^(0 = r 3/2 + 3r 1/2 sinh- 1 £. (32) 

In fig. (18), we plot the ratio of A^^O, 0) for a rectangle to A&^O, 0) for a circle of the 
same area, as a function of the edge length ratio of the rectangle. We see that A fcvR (0, 0) 
exhibits only a mild shape dependence. 

Earlier numerical computations of Ck, R {s) H,|5|JTT|| are all in qualitative agreement with 
our considerations; specifically, the average discrepancy between Ck,R(s) and (Ck,R(s)) is 
always roughly given by (/c 2 /^) -1 / 4 . A detailed comparison is hindered by two issues. 
First, in accord with Berry's original definition, earlier authors usually work with an 
autocorrelation function (7 fejB (s) which, in our notation, is 

From our point of view, Ck,R{s) is a much more complicated object than Ck,R(s); there is no 
simple expression for (C^r(s)), because the relevant functional integral (of Ck,R (s) / Ck,R(0) 
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times the probability distribution P(ipk) over ipk) cannot be done using the simple gaussian 
combinatorics of (flO|) . Of course, by definition (7^(0) = 1, and so essentially what happens 
is that any fluctuation in 6^(0) shows up as a multiplicative enhancement or suppression 
of fluctuations in Cfc,i?( s ) at nonzero s. 

Another problem occurs if an axis of symmetry of the billiard passes through the 
averaging region. Every energy eigenfunction is either symmetric or antisymmetric under 
reflection about such an axis; this can be handled analytically by writing 



1 

71 



Xk{x,y) + Xk(x, -y) 



(34) 



(where we have illustratively assumed an eigenfunction symmetric about the x-axis), and 
treating Xfc( x ) as a gaussian random variable. This approach considerably enhances the 
complexity of the analysis, however; for example, the number of independent terms on 
the right-hand-side of (^) grows from three to forty-eight. The simplest solution is to do 
numerical analysis with averaging regions that do not cross any axes of symmetry, as we 
have done here. 

Finally, we note that all of our results have a straightforward generalization to higher 
dimensions. For a D-dimensional billiard, the autocorrelation function becomes fllll 



(Cfc,fl(s)> 



Jd D p5(p 2 - k 2 )e ip - s 
J d D p 5{p 2 - W) 

2 (D-2)/2 T( n/2) J(D-2)/2(ks) 

[ 1 ' {ks)^ 2 )/ 2 



EE F D {ks) , 

where J v {x) is a Bessel function. The generalizations of 
the properties of Fz>(ks), are 



Ol) and 



(35) 

which follow from 



Afc r(s!, S 2 , 



Id 



D-UJD-l)/D 



F D (k\ Sl - s 2 \) + F D (k\ Sl + s 2 \) 



and 



Afcp(si,s 2 ) 
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D 



k 



D- 



l V {D ~ 

V R 



X)/D 



F D {k Sl )F D {ks 2 ) 



(36) 



(37) 



where Vr is the D-dimensional volume of the spherical averaging region, and 7/5 is a 
numerical factor which we have not computed. 

To conclude, we have performed an analysis of the autocorrelation function Ck } R{s) under 
the assumption that the energy eigenfunction ^(x) behaves like a gaussian random variable, 
in a sense which we have made precise. We find that, for a two-dimensional billiard, Ck t ii(s) 
should have 0{% 1 ^ 2 ) fluctuations about its expected value (Ck t ii(s)) = Jo(ks); scars from 
isolated periodic orbits would give corrections to Cfc^(s) which are also Oi^h 1 ^ 2 ). We have 
given analytic formulae for the root-mean-square amplitude of the expected fluctuations 
in Ck,n(s). We find that a particularly useful object to study is C^^s), which is C^^s) 
averaged over the angle of s. We predict that Ck,R(s) / Jo(ks) is independent of s, a prediction 
which is very well satisfied by the numerical results of Li and Robnik for the Robnik billiard. 
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FIGURE CAPTIONS 

Fig. 1 The autocorrelation function C\r(s) is shown as a solid line. The gray band 
depicts (C ktR (s)) ± A fcjfi (s, s) 1 / 2 , which is the expected root-mean-square range of C ktR (s). 
The Robnik billiard is shown in the inset; the averaging region R is indicated by the filled 
circle, and the direction of the separation vector s is indicated by the direction of the two- 
headed arrow, which has unit length. 

Figs. 2-12 Same as fig. (1). 

Fig. 13 The autocorrelation function Ck,R(s), averaged over the direction of the 
separation vector s, is shown as a solid line. The gray band depicts (Cfc^s)) ± A^(s, s) 1 / 2 , 
which is the expected root-mean-square range of Cfc^s). The Robnik billiard is shown in 
the inset; the averaging region R is indicated by the filled circle. 

Figs. 14-16 Same as fig. (13). 

Fig. 17 [Ck,R.(s) — Ck,R(0)Jo(ks)] + Cfc ; ft(0) for each of the four averaging regions. This 
quantity is predicted to be independent of s. 

Fig. 18 The ratio of A fej #(0, 0) for a rectangle with edge length ratio b/a to A fejfl (0, 0) 
for a circle of the same area. 
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